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Abstract 

In this study we investigate the statistics of two-dimensional stationary turbulence using a Markovian forcing 
scheme, which correlates the forcing process in the current time step to the previous time step according to 
a defined memory coefficient. In addition to the Markovian forcing mechanism, the hyperviscous dissipation 
mechanism for small scales and the Ekman friction type of linear damping mechanism for the large scales 
are included in the model. We examine the effects of various dissipation and forcing parameters on the 
turbulence statistics in both wave space and physical space. Our analysis includes the effects of the effective 
forcing scale, the bandwidth of the forcing, the memory correlation coefficient, and the forcing amplitude, 
along with the large scale friction and small scale dissipation coefficients. Scaling exponents of structure 
functions and energy spectra are calculated, and the role of the parameters associated with the Markovian 
forcing is discussed. We found that the scaling exponents are approximately invariant and show a universal 
behavior for the various forms of forcing schemes used. We found, however, that the final states strongly de- 
pend on the large scale friction mechanism considered. When the large scale friction mechanism is included 
in the model with a high friction coefficient, we demonstrate that the behavior is no longer universal. Our 
analysis also shows that the second-order vorticity structure function has an asymptotic scaling exponent for 
larger dissipation. Additionally, we confirmed that vorticity behaves as a passive scalar when the dissipation 
mechanism becomes less effective. Finally, although turbulence is not believed to have a separation of time 
scales in the dynamics of the velocity field, we conjectured that a separation of time scales exists in the 
dynamics of the energy spectrum. 
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1. Introduction 



The phenomenology of turbulence was described by Richardson [1 and quantified in a scaling theory 
by Kolmogorov [2 . It is believed that the turbulent fiow phenomena are describable through the three- 
dimensional (3D) Navier-Stokes equations Two-dimensional (2D) turbulence, to the first approximation, 
is a reduced dimensional version of 3D turbulence, where the fiow is constrained to two dimensions. In reality, 
2D turbulence is never realized in nature or in the laboratory, both of which have some degree of three- 
dimensionality [4 . Nevertheless, many aspects of idealized 2D turbulence appear to be relevant for physical 
systems in geophysics, astronomy and plasma physics [5 . One of the most important reasons for studying 
two-dimensional turbulence is to improve our understanding of geophysical fiows in the atmosphere and 
ocean [6HT3], in which stratification and rotation suppress vertical motions in the thin layers of fiuid. 

From a theoretical perspective, 2D turbulence is not simply a reduced dimensional version of 3D turbu- 
lence because a completely different phenomenology arises from new conservation laws in two dimensions 
|4]. In fact, two-dimensional turbulence behaves in a profoundly different way from three-dimensional tur- 
bulence due to different energy cascade behavior, and follows the Kraichnan-Batchelor-Leith (KBL) theory 
p!4]4T6) . In three-dimensional turbulence, energy is transferred forward, from large scales to smaller scales, 
via the vortex stretching and tilting mechanism. In two dimensions that mechanism is absent, and it turns 
out that under most forcing and dissipation conditions energy will be transferred from smaller scales to 
larger scales. This is largely because of another quadratic invariant, the potential enstrophy, defined as 
the integral of the square of the potential vorticity. Despite the apparent simplicity of dealing with two 
rather than three spatial dimensions, two-dimensional turbulence is possibly richer in its dynamics due to 
its conservation properties, such as its inverse energy and forward enstrophy cascading mechanisms, which 
three-dimensional turbulence does not possess. 

The physics of two-dimensional turbulence have been elucidated substantially during the past decades by 
theoretical models, intensive numerical investigations, and dedicated soap film experiments [171 [18]. Danilov 
and Gurarie [19 , Kellay and Goldburg [20 , and Tabeling [21 reviewed both theoretical and experimental 
two-dimensional turbulence studies. More recent reviews on two-dimensional turbulence are also provided 
by Clercx and van Heijst [22 and Boffetta and Ecke [4 . Recent studies in two-dimensional turbulence, both 
forced (stationary) turbulence [23H28] and unforced (decaying) turbulence [2'9U32| provide high resolution 
computational confirmation of the KBL theory. The conjecture in KBL theory is that enstrophy, not energy, 
cascades to the small scales, and the energy, on the contrary, cascades to the large scales [33 . Therefore, 
in two-dimensional turbulence there are two inertial ranges, one for the forward cascade of enstrophy and 
one for the inverse cascade of energy. In this dual cascading phenomenon, the relative locations of the 
Kolmogorov and Kraichnan scalings in 2D turbulence energy spectra depend on the forcing scale [34 . 
According to Kolmogorov theory, in the energy cascade range, the only parameters of practical importance 
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would be energy injection rate e and wave number k. Dimensional reasoning states that the energy density 
(i.e., energy spectrum) is E(k) ^ e^/^/c~^/^. In KBL theory, a similar argument gives that E{k) ~ rf^^k~^ 
in the forward enstrophy cascade range; where r] is the enstrophy injection rate (i.e., energy injection rate 
is related to the enstrophy injection rate by e = rj/k'^^ in which /c/ is the energy injection scale). These 
predictions for dual cascade given by KBL theory are illustrated in Figure [l] 

inverse energy 




' small scale 

I dissipatio n 

Waven umber 

Figure 1: A qualitative picture of the energy spectrum which shows double cascade scalings of stationary two-dimensional 
turbulence. 

The dual cascading conjecture, along with the dissipation mechanisms, challenges the universality of 
the scalings in stationary two-dimensional turbulence. In general, the dissipation mechanism contains an 
energy sink on large-scales and an enstrophy sink on small-scales [lOl HH [35H39] . The presence of these 
two sinks is necessary to reach a stationary regime if two different fluxes are assumed to flow in opposite 
directions from the forcing zone as hypothesized in two-dimensional turbulence. In other words, due to 
inverse energy cascading in two-dimensional turbulence, the stationary final state can only be obtained if 
there exist a dissipation mechanism in large scales. Therefore, large-scale dissipation mechanisms have been 
routinely employed in numerical simulations of stationary two-dimensional turbulence to absorb energy at 
large scales. Blackbourn and Tran [40] have recently studied the effects of friction on forced two-dimensional 
Navier- Stokes turbulence. 

Two-dimensional turbulence models are useful for understanding large scale motions of forced-dissipative 
systems occurring in the atmosphere and oceans. The frictional and forcing effects discussed in this study 
are crucial to understanding the dynamics of these systems. The statistics of such a forced-dissipative two- 
dimensional turbulence system are investigated here for a wide range of physical considerations. One of our 
major goals is to analyze the universality of the scaling laws for these statistics in both wave space and 
physical space by considering the energy or enstrophy spectra and velocity or vorticity structure functions. 
In order to minimize the effects of dissipation at intermediate scales of the spectrum, following many two- 
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dimensional turbulence studies (e.g., [24 l I35H39] ) . we use high-order Laplacian for viscosity (sometimes 
called hyperviscosity) which separates sharply the inertial and dissipative ranges. To model the forcing 
and energy injection mechanism, we use a Markovian process, which correlates the forcing structure in the 
current time step to the previous one according to a defined memory coefficient. The effects of the amplitude 
of the forcing, the effective forcing scale, and the bandwidth of the forcing are considered here in addition 
to the memory coefficient parameter that measures the stochastic behavior of the forcing mechanism in 
the system. We integrate the hyperviscous Navier-Stokes equations with a pseudospectral method using 
the third-order Runge-Kutta scheme for nonlinear terms and the second-order Crank-Nicolson scheme for 
the linear dissipation terms. We should also note that our hyperviscous model reduces to the classical 
Navier-Stokes turbulence model in the limiting case. 

The paper is organized as follows: the mathematical formulation of the forced-dissipative two-dimensional 
turbulence is given in Section [2] The numerical methods are brieffy presented in Section [3) The results for 
two-dimensional isotropic homogeneous stationary turbulence are provided in Section [4] This section also 
explores the effects of large scale friction and small scale dissipation mechanisms, and systematically analyzes 
the relevant parameters of the Markovian forcing scheme. The turbulence statistics in the stationary regime 
and their scaling exponents are shown in this section as well. Final conclusions and comments are drawn in 
Section [5l 



2. Mathematical model 

The governing equation for two-dimensional incompressible ffows can be written in its vorticity formu- 
lation in the following form \19^ 

%.^^D.F ,1, 

where uj is the vorticity which is a scalar quantity in two-dimensional ffows. Here, J is the nonlinear Jacobian 
which symbolizes nonlinear interactions, D represents the dissipation mechanism, and F is for forcing. The 
Jacobian term in Eq. ([T]) is deffned as 

j = „.Vu;=^^-^^ (2) 

dy dx dx dy 

using the deffnition of the velocity stream function 

u=^' v = -^ (3) 
dy ' dx 

where u and v are components of the two-dimensional velocity vector ffeld u. The dissipation term in 
classical Navier-Stokes equation is given as I) = ^;V^a;, where is the Laplacian. In this study, we use a 
generalized form of the dissipation D = vV'^'^uo in place of the regular viscosity dissipation. The purpose of 
using high order Laplacian, which is called as hyperviscosity, is to eliminate as much as possible the effects 
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of viscosity at the intermediate scales, thus extending turbulence inertial scales. On the other hand, in order 
to sink energy at the large scales we use Ekman type friction to be able reach a statistically steady state. 
In the case of periodic boundary conditions the applied forcing mechanism would result in an unbounded 
growth of the total kinetic energy if we would not include a large scale friction mechanism. The Ekman 
friction terminology is usually used in the contexts of the rotating flows. We can also interpret this large 
scale friction term as Rayleigh friction for stratified flows, or Hartman friction for magnetohydrodynamic 
flows. Finally, the dissipation mechanism in our study is modeled by the following generalized form 

D = {-If+^vW^Puj - Xcu (4) 

in which it reduces to the classical Ekman-Navier- Stokes equation for p = 1. The kinematic relationship 
between vorticity and stream function according is given as 

^ + ^ = _^ (5) 

From a computational point of view, this formulation has several advantages over the primitive variable 
formulation. It eliminates pressure from the Navier-Stokes equations and hence has no corresponding odd- 
even decoupling between the pressure and velocity components, as well as projection inaccuracies usually 
observed in fractional step approaches [41 . The vorticity-stream function formulation also automatically 
satisfies the divergence-free condition and allows one to reduce the number of equations to be solved. 



3. Numerical methods 

Fourier series expansion based methods are often used for solving problems with periodic boundary 
conditions. One of the most accurate methods for solving the Navier-Stokes equations in periodic domains 
is the pseudospectral method, which exploits fast Fourier transform (FFT) algorithms, resulting in spectral 
accuracy g2j [43] . By transforming Eq. ([T]) to Fourier space the governing equation becomes 

^ + Jk = l^k + i^k (6) 

where a hat over the variable represents the Fourier coefficients of the corresponding variable in the wave 
space k = {k^^ ky). The relationship between vorticity and stream function in the Fourier space becomes 

{-kl - kl)i,^ = -Co^. (7) 

The nonlinear Jacobian in Fourier space is 

Jk = {iky2pi^) o {ikxUJ\^) - (ikx'ipi^) o {ikyLU\^) (8) 

where i is the complex unit number (i.e., f = —1). The convolution sum in the nonlinear Jacobian term 
is computed in the spatial domain using the convolution theorem. In the pseudospectral method, the 
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convolution sum of these nonlinear terms is actually computed in the physical domain and Fast Fourier 
transforms are used to go back and forth between Fourier wave space and physical space. In the present 
study, we use the standard Fourier- Galerkin pseudospectral method in a periodic square box of a length 27r 
j44| [45] . The dissipation terms in the Fourier space becomes 

k 



[vk^P + A] 



kd 



(9) 



U2p 



where k = ^k"^ -\- k^. We redefine the small scale viscosity coefficient by using the expression of u ■ 
where k^ is the effective dissipation wave number, which we set in our study as k^ = 0.96A/'/2, where 
N'^ = 512^ is the resolution of the problem. Note that the power p = 1 corresponds the regular constant 
property Navier-Stokes equations for Reynolds number of Re = k'^/u. In our model, we use a Markovian 
forcing scheme, which correlates the forcing in the current time step to the previous one according to a 
deffned memory coefficient. The forcing is localized within narrow spectral range {kf — cr^kf -\- a) in the 
vicinity of the forcing wave number kf. The bandwidth of the effective forcing is determined by the a 
variable. The Markovian forcing process at the current time step becomes 



f,il-py/^e<+pfr' 



(10) 



where /o is the forcing amplitude, p is the memory correlation coefficient, and the ( is the uniformly 
distributed phase on interval [0,27r]. In this study, the semi-discrete vorticity transport equation in Fourier 
space, Eq. ([6|, is solved by a combination of the third-order Runge-Kutta and second-order Crank-Nicolson 
schemes in a periodic square. The dissipation terms are treated by an implicit Crank-Nicolson scheme and 
nonlinear Jacobian and forcing terms are treated by the explicit Runge-Kutta scheme. Therefore, starting 
with the value of the Fourier coefficients of vorticity, at the current time step, the time marching 
algorithm for computing the vorticity at the next time step, consists of the following three substeps: 



^(2) 



\n+l 



(2) 



■72(-J^ + Fk") 
■73(-4'^+i^k")- 
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where the coefficients are 
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in which At is the time step. 
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4. Results 



In this section we present numerical results for homogeneous isotropic stationary turbulence for various 
physical parameters which determine the forcing and dissipation mechanisms. The spatial resolution is fixed 
N'^ = 512^ Fourier modes, the time step size is At = 0.005, and the solution is advanced forward in time 
in Fourier space. For all the simulations we start from a rest state, integrate the model until a statistically 
steady state is obtained, and continue for enough time to compute turbulence statistics. The boundary 
conditions are 27r-periodic in both directions. The domain-integrated kinetic energy is tracked for each 
simulation to measure the energy level of the system, which is quantified by the following integral 

in which we can estimate when the system reaches quasi-stationary regime. Averaging over time is not 
possible unless a statistically steady state is established by a large scale friction mechanism. The mean 
values are obtained by ensemble averaging the data between time t = 50 and t = 100 throughout the study. 

In order to examine the characteristics of two-dimensional stationary turbulence, we first define two 
statistical measures; one is the energy spectrum in wave space, and the other is the structure function in 
physical space. The energy spectrum is defined as 

E{k,t) = ^k^\i>{k,t)f (14) 

and the angle averaged energy spectrum is 

E{k,t)= J2 ^('^'^)- (1^) 

/c<|k|</c+l 

In the forward enstrophy cascade (i.e., in the inertial scale between the forcing scale and the small scale 
dissipation) , it is known from the KBL theory that the energy spectrum in the inertial range approaches the 
classical scaling in the inviscid limit. On the other hand, in the inverse energy cascade (i.e., the inertial 
scale which includes the scales larger than the forcing scale), the theory predicts that the energy spectrum 
scales by A:~^/^, since the scaling arguments leading to this Kolmogorov scaling nowhere assume a specified 
direction of the cascade. 

The statistics of two-dimensional turbulent flow can be further investigated considering powers of velocity 
or vorticity diff'erences in the physical space. A commonly used statistical quantity in two-dimensional 
turbulence is the second-order vorticity structure function which is deflned as 

(M^)2) = (mx + r)-c.(x)|2) (16) 

with r = |r| being the spatial separation. Assuming that the system is homogeneous and isotropic, the 
structure functions depend on r only. The classical k~^^^ scaling of energy spectrum in the inverse energy 
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Figure 2: The effects of the small scale dissipation coefficient on the statistics (A = 0.05, p = 8 and kf = 5); (a) time series of 
total energy, (b) angle averaged energy spectra, and (c) second-order vorticity structure functions. 
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Figure 3: Instantaneous vorticity fields at time t = 100 for varying the small scale dissipation coefficient using the forcing scale 
kf = 5 and the order of Laplacian p = 8. 



cascade relates to a r^/^ scaling of second-order velocity structure function. In the enstrophy cascade range, 
where Kraichnan scaling appears, the corresponding scaling is equal to the upper bound. Hence, 
steeper energy spectra always results in scaling [46 . Vorticity transport equation in two-dimension is 
formally identical to passive scalar transport equation that describes the transport of a scalar quantity 
6 passively advected by velocity field ^47H49] . In the context of passive scalar, according to the theory 
the second-order vorticity structure function in the inviscid limit is proportional with (Jcj(r)^) ~ 
for a given energy spectrum E{k) ~ k~^~^^ for < (i < 2. 

The scaling exponents of turbulence statistics that deviate from the values predicted by the theory have 
been observed in the literature |^ ITOl [T9l l38| WJ\ [5TH53] . The universality of the statistics of two-dimensional 
turbulence in both wave space and physical space have also been questioned [54 . In the following analysis, 
statistical properties of the forced two-dimensional turbulence are systematically investigated by numerical 
simulations for different physical parameters. Our primary goal here is to elucidate the effects of the various 
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Figure 4: The effects of the order of hyperviscosity on the statistics (A = 0.05, = 10 and kf = 5); (a) time series of total 
energy, (b) angle averaged energy spectra, and (c) second-order vorticity structure functions. 

forcing and dissipation mechanisms on the turbulence statistics. 

4.1. Effects of small scale dissipation mechanism 

First, we analyze the effects of small scales dissipation mechanism on turbulence statistics. Figure [2] 
shows the statistics for different values of small scale dissipation coefficient u. The order of the Laplacian, 
the large scale dissipation coefficient, and all the coefficients for forcing scheme remain fixed {p = 8, A = 0.05, 
/c/ = 5, cr = 2, and p = 0.5). As we can see from the figure, an increase in the dissipation coefficient results 
in a slightly smaller inertial range in the energy spectrum; however, the scalings of the spectra remain the 
same. This deviation can also be seen from the comparison plot of the vorticity structure function for smaller 
scales. The levels of energy of the systems in the stationary regime are close to each other. This comparison 
clearly shows that the small scale dissipation mechanism has slight effect on the statistics, due to the sharp 
effect of the high order Laplacian in the dissipation mechanism. The instantaneous vorticity fields are also 
plotted in Figure [3j showing that there are no significant effects on flow field structures as well. 

Second, the effects of the order of the Laplacian in the small scale dissipation mechanism are investigated 
by using the same coefficients given above for other mechanisms (A = 0.05, kf = 5, a = 2, and p = 0.5). The 
case for the order of Laplacian p = 1 represents regular quadratic dissipation for Navier-Stokes equations, 
and all the other cases of p > 1 represent the hyperviscous Navier-Stokes equations. As shown in Figure [4j 
increasing the hyperviscosity order p extends the inertial range and results in the small scale dissipation 
becomes effective in larger wave numbers. Its effect on the flow field is also illustrated in Figure [5| showing 
that the p = 1 case has a less amount of filamentation compared to the hyperviscous cases. The comparison 
of vorticity structure functions in Figure |4] also demonstrates that the scaling for the case of p = 1 is 
{Suj{r)'^) ~ for smaller separation distances. It can also be seen that the vorticity structure functions 




(a) p = 1 (h) p = 2 (c) p = 4 (d) p = 6 

Figure 5: Instantaneous vorticity fields at time t = 100 for various order of the hyperviscosity using the forcing scale of kf =5 
and the small scale dissipation coefficient of u = 10. 




(a) (b) (c) 

Figure 6: The effects of the small scale dissipation coefficient on the statistics without hyperviscosity (A = 0.01, p = 1, p = 0, 
a = 2, fo = 0.1 and kf = 5); (a) time series of total energy, (b) angle averaged energy spectra, and (c) second-order vorticity 
structure functions. 

flatten for large separation distances as predicted by the KBL theory of two-dimensional turbulence in the 
inviscid limit. 

Next, we investigate the Reynolds number dependence of the statistics of stationary turbulence using the 
the order of Laplacian p = 1 that represents the classical quadratic dissipation mechanism in Navier-Stokes 
equation without using hyperviscosity. As demonstrated in Figure [6| the angle averaged energy spectrum 
asymptotically reaches the scaling in the inertial range as u decreases (i.e., Re increases). We find that 
the Reynolds number dependency is more stringent if we look at the turbulence statistics in wave space 
using the angle averaged energy spectrum. The structure functions are proportional to for the smaller 
separations and fiatten for higher separations. The corresponding instantaneous vorticity fields at time 
t = 100 are compared in Figure [7| for the same set of Reynolds numbers. As we can see from Figure [7| 
the amount of filament at ion increases for higher Reynolds numbers. Due to the smaller convection in lower 
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(a) 1/ = 5 {Re = 12080) (b) = 10 {Re = 6040) (c) u = 2h {Re = 2416) (d) u = 100 {Re = 604) 

Figure 7: Instantaneous vorticity fields at time t = 100 for varying Reynolds number using the forcing scale kf = 5 and the 
order of Laplacian p = 1. 




(a) (b) (c) 

Figure 8: The effects of the large scale dissipation coefficient on the statistics {kf = 20, /o =0.1, a = 2, p = 0.0, v = 1000 and 
p = 8)] (a) time series of total energy, (b) angle averaged energy spectra, and (c) second-order vorticity structure functions. 

Reynolds numbers, the interaction between two vortices is not as strong as that of the computations with 
higher Reynolds numbers. 

4-2. Effects of large scale dissipation mechanism 

In order to investigate the effects of large scale dissipation mechanism in the dynamics of two-dimensional 
homogeneous turbulence, we perform a similar analysis by varying the large scale friction coefficient, while 
the other parameters associated with small scale friction and forcing mechanisms are held constant in the 
model. Figure [S] shows the computed statistics for the effective forcing scale kf = 20 applied in all cases. The 
statistics on the evolution of the total energy clearly demonstrate that the main mechanism that determines 
the time scale for stationary regime is the large-scale friction mechanism. Due to the inverse energy cascading 
in two-dimensional turbulence, a statistically steady state is established by damping mechanism in low- 
wavenumbers. Increasing the large scale damping coefficient A provides an earlier statistically steady state 
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(a) A = 0.005 



(b) A = 0.01 



(c) A = 0.02 



(d) A = 0.03 



Figure 9: Instantaneous vorticity fields at time t = 100 for varying large scale friction coefficient A using the forcing scale is 
kf = 20 and the small scale dissipation coefficients are u = 1000 and p = 8. 
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Figure 10: The effects of the large scale dissipation coefficient on the statistics (kf = 10, fo = 0.1, a = 3, p = 0.0, u = 1000 and 
p = 8); (a) time series of total energy, (b) angle averaged energy spectra, and (c) second-order vorticity structure functions. 



with a less amount of energy. It is interesting to note that the structure of the large scale spectrum 
also depends considerably on the amount of the large scale friction coefficient in which Kolmogorov scaling 
appears for smaller value of the A. It can also be seen from Figure [S] that energy spectra in forward enstrophy 
cascade range appear steeper than with increasing large scale damping coefficients. We found that the 
vorticity structure functions are independent of the separation length for larger r, as indicated by KBL 
theory in the inviscid limit, and change gradually from to r^/^ for smaller separations. Thus, we conclude 
that the scaling exponents of the structure function are considerably influenced by the large scale damping 
mechanism for the small separation r. The instantaneous vorticity fields compared in Figure [9] show another 
interesting observation that the flow pattern has more vorticity filaments for smaller large scale friction 
coefficients, whereas the amount of vortical structures seems equivalent in each case due to the use of the 
same forcing scale. 



Figure 10 shows statistics for /c/ = 10. Similar to the previous comparisons with kf = 20, the dual 
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(a) A = 0.005 



(b) A = 0.01 



(c) A = 0.025 



(d) A = 0.05 



Figure 11: Instantaneous vorticity fields at time t = 100 for varying the large scale friction coefficient A using the forcing scale 
kf = 10 and the small scale dissipation coefficients of u = 1000 and = 8. 
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Figure 12: The effects of the large scale dissipation coefficient on the statistics (kf = 5, fo = 0.1, a = 2, p = 0.5, u = 1000 and 
p = 8)] (a) time series of total energy, (b) angle averaged energy spectra, and (c) second-order vorticity structure functions. 



cascading phenomenon, Kraichnan scaling for the scales smaller than the energy injection scale, and Kol- 
mogorov scaling for the scales greater than the forcing scale, appears in statistics for energy spectra and 
structure functions. The lines k~^^^ and are shown in the plot for energy spectra and lines r^/^ and 

are included in the plot for structure functions for comparison. The corresponding flow patterns are also 
shown in Figure [TT] demonstrating a clear comparison for the effect of the large scale friction coefficient into 
the flow structure. For larger friction coefficient, flow tends to be concentrated in the centers of correspond- 
ing vortices, showing less interactions among them. The grid-like pattern of vortices emerges, corresponding 
to the accumulation of energy near forcing scales; this pattern is also evident in energy spectra. 

In order to further investigate the effects of large scale friction mechanism, we next perform a similar 
analysis for a larger energy injection scale kf = 5. The results are summarized in Figure 12, showing 



evolution of total kinetic energy, mean energy spectra, and mean second-order vorticity structure functions 
for a series of runs with varying large scale friction coefficient. For A = 0, the flow approaches two- 
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(a) A = 0.0 (b) A = 0.01 (c) A = 0.025 (d) A = 0.1 

Figure 13: Instantaneous vorticity fields at time t = 100 for varying the large scale friction coefficient A using the forcing scale 
kf = 5 and the small scale dissipation coefficients of u = 1000 and p = 8. 

dimensional homogenous turbulence with a forward cascade of enstrophy and an inverse cascade of energy. 
After an initial period of nonlinear adjustment, the total energy grows linearly in time, dE/dt ~ t. Some 
part of the energy input injected by forcing mechanism at forcing scale kf = 5 is transferred to scales 
smaller than the forcing scale, and dissipated by viscosity. In this finite system for A = 0, however, there 
is no statistically steady state due to inverse energy cascade in which energy is continuously transferring 



from the forcing scale to the larger scales. Figure 12 also shows that as the large scale friction coefficient is 



decreased, the growth of energy is decreased, until at a value A ~ 0.025, the growth is completely suppressed 



by the large scale damping, and stationary turbulence is obtained. Comparing energy spectra in Figure 12 
Kraichnan scaling is obtained for A ^ 0, though the spectra for larger A appear steeper than k~^. It is also 
important to note that the tails of corresponding second-order vorticity structure functions scale as for 
larger A. Therefore, we confirm that the scaling for the structure functions shows an asymptotical limit for 
ffows having energy spectrum is steeper than k~^. Based on these findings, we conjecture that looking for 
turbulence statistics only in terms of structure functions does not guarantee the proper scaling exponents 
in energy cascading predictions. Figure [12] also shows that the second-order vorticity structure function 



fiattens for larger separation distances r, as predicted by the KBL theory in the inviscid limit. Figure 13 
illustrates the instantaneous vorticity field at time t = 100 for a set of values for A. Similar to our previous 
outcomes, we observe more layered fiow patterns for smaller A, while energy is concentrated in the vortical 
structures associated with the forcing scale for larger A. 



Figure [M] shows energy spectra at several times to illustrate the spectral energy density evolution for a 
series of cases with two values of A and kf. The time interval between two adjacent lines is 5t = 2. The 
damping effect of the large scale friction mechanism in the development of the energy cascade is observable 
for the kf = 5 cases in the early time evolution. We also observe that a separation of time scales exists in the 
dynamics of the energy spectrum. This can be seen more clearly for the cases in which kf = 15. Once the 
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(a) kf =5, X = 0.0 



(b) kf =5, X = 0.01 




(c) kf = 15, A = 0.0 (d) kf = 15, X = 0.01 

Figure 14: Evolution of the angle averaged energy spectrum (/o =0.1, a = 2, p = 0.5, u = 1000 and p = 8) for varying the 
large scale friction coefficient A and the effective forcing scale kf. 

forcing of the initially at-rest fluid begins, the entire spectrum immediately fills out with noise and a strong 
forcing peak appears centered around a wave number commensurate with the forcing length scale. Certain 
wave number ranges in a typical energy spectrum in a developing turbulent fiow are populated immediately, 
while certain characteristics, the k~^^^ scaling in the wave numbers smaller than the forcing scale, for 
example, evolve slowly over time. The spectrum fills out slowly until, after a long time, the final statistically 
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steady state Kolmogorov spectrum appears if there exist even a small amount of large scale damping. This 
may be contrasted with the time evolution of a three-dimensional homogeneous isotropic turbulent flow, 
which develops its structure immediately, and with time transfers this structure to progressively smaller 
length scales. 

4-3. Effects of forcing mechanism 

Here, we systematically analyze the effects of the parameters associated with Markovian forcing scheme. 
The underlying settings for the small scale dissipation mechanism are same for all the cases (i.e., u = 1000 
and p = S). The use of hyperviscosity maintains a constant flux of enstrophy in a wider interval kf > k > 
by eliminating as much as possible the effects of viscosity at intermediate scales, thus extending the inertial 
ranges. First, we compute the statistics by varying the corresponding scale for Markovian forcing using a 



large scale friction coefficient of A = 0.05. Figure 15 shows the statistics in terms of evolution of the total 



energy, mean energy spectrum, and mean second-order vorticity structure function for different forcing scales 
with the same memory correlation coefficient, forcing amplitude, and forcing bandwidth. The comparison 
of time series clearly demonstrates that all computations reach quasistationary regime at the same time but 
having a different level of energy. Increasing the forcing scale result in a decrease in quasistationary energy 
level. It is also shown that the variability increases by decreasing forcing scale. On the other hand, the 
scaling exponents in energy spectra and structure functions are invariant of the forcing scale. Instantaneous 



vorticity fields at time t = 100 for varying the forcing scale kf are also illustrated in Figure 16 showing the 
similar pattern of mixing and filamentation. 

Similarly, we compute the statistics by varying the effective scale for Markovian forcing using a smaller 



large scale friction coefficient of A = 0.01. The statistics and flow field patterns are shown in Figure \17\ and 
Figure [Tsj respectively. Due to the reduction in damping coefficient we observe more interaction between 

16 



(a) kf = 3 



(b) kf = 5 



(c) kf = 10 



(d) kf = 15 



Figure 16: Instantaneous vorticity fields at time t = 100 for varying the forcing scale kf using the large scale dissipation 
coefficient of A = 0.05. 
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k = 10 
k = 15 
k = 20 



(b) 



k = 5 
k = 10 

k = ^5 

k=20 
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Figure 17: 
energy, (b) 



The effects of the forcing scale on the statistics (A = 0.01, a = 2^ p = 0.0, and /o 
angle averaged energy spectra, and (c) second-order vorticity structure functions. 



0.1): 



(c) 

(a) time series of total 



vortical structures for all the cases with varying kf. In the inertial range, especially for /cj = 5, it is also 
interesting to see that energy spectrum scales as while structure functions scales as r^/^. The main 
reason for energy spectrum appearing steeper than for increasing kf \s that the relative importance of 
the large scale damping coefficient in larger kf cases is greater. This can also be seen from the energy levels 
shown in Figure [T7[ 

Next, we analyze the effects of the memory coefficient p in the Markovian forcing scheme. As shown in 
Figure [T9j its main effect on the statistics is the translation of the level of energy. There is no significant 
effect into the scaling exponents of the energy spectra and structure functions. The vorticity fields are also 



illustrated in Figure 19 showing that the amplitude of vorticity increases with increasing memory coefficient 
p. In this study, we focus on the random forcing mechanism controlled by this coefficient that measures the 
stochastic process in the system for which the forcing is purely random (i.e., p = 0). It will be interesting 
to further investigate the effects of deterministic forcing in the context of generation of coherent vortices. 
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(a) kf = 5 



(b) kf = 10 



(c) kf = 15 



(d) kf = 20 



Figure 18: Instantaneous vorticity fields at time t = 100 for varying the forcing scale kf using the large scale dissipation 
coefficient of A = 0.01. 
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p = 0.25 
p = 0.5 
p = 0.75 
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(c) 



Figure 19: The effects of the memory coefficient on the statistics (A = 0.05, kf = 5, a = 2, u = 1000 and p = 8); (a) time series 
of total energy, (b) angle averaged energy spectra, and (c) second-order vorticity structure functions. 



However, we particulary use a random forcing mechanism to concentrate on the statistical deviations from 
the theoretical scalings which are solely due to viscous effects and eliminate possible situations that coherent 
vortices might play a role. 



Figures 21 and 22 examine the effects of forcing bandwidth a for the forcing scale kf = 10. The effective 



energy injection occurs in the scales within the interval {kf — a, kf -\-a). As shown in Figure 21 , the level of 
energy of the system in statistically steady state increases with increasing the a since we inject more energy 
due to enhanced forcing zone. In the forward energy range we observe the same scaling for energy spectrum. 
For cr = 8, it is also shown that the energy spectrum within the forcing zone scales as E{k) ~ k~^^^. 
Comparing the structure functions, we illustrate that {5uj{r)'^) wiggles more rapidly for smaller values of 
<j. It is also interesting to see that the flow pattern shows more vortical filaments for increasing forcing 
bandwidth. 



Finally, we show the effects of forcing amplitude /o in Figures 23 and 24 The linear translation of 

18 




(a) p = 0.0 (b) p = 0.25 (c) p = 0.5 (d) p = 0.9 

Figure 20: Instantaneous vorticity fields at time t = 100 for varying the memory coefficient using the large scale dissipation 
coefficient of A = 0.05 and the small scale dissipation coefficients oi u = 1000 and p = S. 




(a) (b) (c) 

Figure 21: The effects of the forcing bandwidth on the statistics (A = 0.05, kf = l^,p = 0.5, u = 1000 and p = 8); (a) time 
series of total energy, (b) angle averaged energy spectra, and (c) second-order vorticity structure functions. 
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(a) 



(b) 



(c) 



Figure 23: The effects of the forcing amphtude on the statistics (A = 0.05, kf = 5, p = 0.5, u = 1000 and p ■ 
series of total energy, (b) angle averaged energy spectra, and (c) second-order vorticity structure functions. 



(a) time 




(a) /o = 0.05 



(b) /o = 0.1 



(c) /o = 0.2 



(d) /o = 0.3 



Figure 24: Instantaneous vorticity fields at time t = 100 for varying the forcing amplitude using the large scale dissipation 
coefficient of A = 0.05 and the small scale dissipation coefficients of u = 1000 and p = 8. 



system's energy level in the statistically steady state can be clearly seen from Figure |23) Here /o represents 
the energy injection rate of the external forcing. Therefore, increasing /o also results in an increase for the 



amplitude of vorticity which can be seen Figure 24 However, the shapes of corresponding energy spectra 
and structure functions look like similar and there is no significant differences for these statistical quantities. 



5. Summary and conclusions 

Numerical simulations of forced two-dimensional turbulence have been performed by solving the Ekman- 
Navier-Stokes equations using the Fourier- Galerkin pseudospectral method. Large scale friction, high order 
small scale dissipation, and Markovian forcing mechanisms have been included in the model. The formulation 
presented here can be reduced to the classical Navier-Stokes model as a special case by setting the Ekman 
friction coefficient equal to zero and the order of the viscosity equal to unity. Including the large scale 
friction mechanism in the model is crucial to be able to obtain a stationary turbulent ffow regime in the case 
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of periodic boundary conditions due to the inverse energy cascading phenomenon. Otherwise, the apphed 
forcing mechanism would result in the unbounded growth of the total kinetic energy. The objective of 
this study was to determine the turbulence statistics for the long-time integration of stationary turbulence 
simulations, and to investigate the scaling exponents for a large range of physical settings. Specifically, we 
computed the statistics in terms of the evolution of total kinetic energy, angle averaged energy spectrum, 
and second-order vorticity structure function. 

We first concentrated on the effects of the small scale dissipation mechanism using hyperviscosity, which 
turns on much more abruptly than the gradual increase of normal viscosity at the small scales. We also stud- 
ied the Reynolds number dependence when using classical viscosity with the order of Laplacian p = 1. We 
showed that the predicted energy spectrum asymptotically converged to the theoretical scaling as the 
Reynolds number increased, which is predicted by the KBL theory for forward cascading two-dimensional 
turbulence. We also showed that hyperviscosity effectively eliminates the effects of viscosity at the interme- 
diate scales, thus extending the turbulence inertial range. We demonstrated the shape of the second-order 
vorticity structure function that is proportional to for the smallest separations and flattens when ap- 
proaching the forcing length scale. 

Next, we studied the effects of the large-scale friction mechanism within the turbulence statistics. The 
linear damping mechanism is utilized in the model, with the friction coefficient A varying from zero to 
higher values. These computations revealed that the existence of a statistically steady state is maintained 
by this mechanism and significantly affects the turbulence statistics. We confirmed the classical dual cascade 
picture of two-dimensional stationary turbulence for A ^ 0, which indicates that a direct enstrophy cascade 
is developed from the forcing range to the dissipation range, and an inverse energy cascade is developed 
for the scales greater than the effective forcing scale. We found, however, that the large scale damping 
mechanism is a major source of deviations from the classical scaling in the forward, and k~^^^ scaling in 
the inverse, cascade ranges. We also showed that the tails of corresponding second-order vorticity structure 
functions scale as for large A, and reduce to the r^/^ scaling as A vanishes. 

Finally, we performed a detailed study by systematically varying the parameters associated with Marko- 
vian forcing mechanism such as the energy injection scale, the forcing bandwidth, the forcing amplitude, 
and memory correlation coefficient. We found that these parameters exhibit no significant difference on 
turbulence statistics, except their translational effects on the total levels of energy in the statistically steady 
state. We demonstrate that the flow patterns show more vortical filaments for increased forcing bandwidth. 
We also showed that there is a separation of time scales in the dynamics of the energy spectrum, such 
that energy is transferred quickly in the forward enstrophy cascade range, while the shape of the spectrum 
fills out slowly, until, after a long time, the final statistically steady state Kolmogorov spectrum appears in 
the inverse energy cascade range. This separation of time scales in the energy spectrum might be a useful 
starting point for developing multiscale computational algorithms for turbulence research, a topic we intend 
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to investigate further in a future study. 
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